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. Abstract 

O . The aim of this Short Note is to highlight that the generalized grey 

' Brownian motion (ggBm) is an anomalous diffusion process driven by 

ly-j . a fractional integral equation in the sense of Erdelyi-Kober, and for 

this reason here it is proposed to call such family of diffusive processes 
as Erdelyi-Kober fractional diffusion. The ggBm is a parametric class 
of stochastic processes that provides models for both fast and slow 
anomalous diffusion. This class is made up of self-similar processes 
' with stationary increments and it depends on two real parameters: < 

& . a < 2 and < (3 < 1. It includes the fractional Brownian motion when 

< a < 2 and (3 — 1, the time- fractional diffusion stochastic processes 
when < a = (3 < 1, and the standard Brownian motion when a = 
/3 = 1. In the ggBm framework, the Mainardi function emerges as a 
^ . natural generalization of the Gaussian distribution recovering the same 

' key role of the Gaussian density for the standard and the fractional 

, Brownian motion. 
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Diffusion in disordered media is commonly named anomalous, to mark a 
difference with the classical diffusion processes where the probability density 
function pdf to find a particle in the place x at the time t is Normal, i.e., 
Gaussian, and for this it is called Normal diffusion, or Gaussian diffusion. 
Moreover, beside the non-Gaussian particle pdf, the anomalously is embod- 
ied even by the non-linear growing in time of the variance of the particle 
spreading [8]. The process of anomalous diffusion is referred to as fast dif- 
fusion, when the variance of the particle spreading grows with a power law 
with exponent greater than 1, and it is referred to as slow diffusion, when 
such exponent is lower than 1. 
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In the case of anomalous diffusion, the classical local flux-gradient re- 
lationship does not hold and it is necessary to determine a non-local re- 
lationship. It is well-known that a useful mathematical tool for physical 
investigation and description of non-local and anomalous diffusion is Frac- 
tional Calculus, which is that branch of mathematical analysis dealing with 
pseudo-differential operators interpreted as integrals and derivatives of non- 
integer order [7, 26J. 

Non-locality can be designated in time (time-fractional diffusion) or in 
space (space- fractional diffusion), as well as both in space and time (space- 
time fractional diffusion equation) [15J . Generally, when the fractional differ- 
entiation is considered for the time, then the fractional derivative operator 
is assumed to be in the Caputo or in the Riemann-Liouville sense, when 
the fractional differentiation is considered for the space, then the fractional 
derivative operator is assumed to be in the Riesz-Feller sense. 

Recently, the extension of fractional differential equations to distributed- 
order fractional differential equations has permitted to describe also pro- 
cesses whose scaling law changes in time, see e.g. [T71 l28l [35] . 

Furthermore, under the physical point of view, when there is no sep- 
aration of timescale between the microscopic and the macroscopic level of 
the process the randomness of the microscopic level is transmitted to the 
macroscopic level and the correct description of the macroscopic dynamics 
has to be in terms of the Fractional Calculus [6j. Moreover, fractional inte- 
gro/differential equations are related to phenomena with fractal properties 

A fractional differential approach has been successfully used for mod- 
elling in several different disciplines as for example statistical physics [19] . 
neuroscience [12], economy [29] . control theory [37] and combustion science 
|24j . Further applications of the fractional approach are recently introduced 
and discussed by J. A. Tenreiro Machado [5B] . 

Normal diffusion, or Gaussian diffusion, is a Markovian stochastic pro- 
cess driven by the classical parabolic equation 



dP d 2 P 



dt dx 2 



xen, ten£, (l) 



with initial condition P(x,0) = Pq(x). The fundamental solution of ([I]), 
which is named also Green function, and corresponding to the case with 
initial condition P(x,0) = Pq(x) = 5(x), is the Gaussian density 

/+oo 
x 2 f(x, t) dx = 2 1. 
-oo 

The density function P(x,t) with general initial condition P(x,0) = Pq(x) 
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is related to the fundamental solution f(x,t) by the following convolution 
integral 

/+oo 
/(£,i)Po(z-£R. (3) 
-oo 

In order to generalize the classical Markovian setting to Non-Markovian 
cases, the following integral equation has been introduced by Mura, Taqqu 
and Mainardi [23] 

P(x,t) = P (x) + ^^K[g{t)-g{s) f P d { ^ s) ds , (4) 

where K(t) is a memory kernel and g{t), with g(0) = 0, is a smooth and 
increasing function describing a time stretching. The Green function of dU 
Q(x,t), which is the marginal one-point one-time pdf of the non-Markovian 
diffusion process, turns out to be 

g{x,t) = / f(x,r)h(r,g(t))dT, (5) 
Jo 

where f(x,t) is the Gaussian density ([2]) that is the fundamental solution of 
the Markovian diffusion process, i.e., K(t) = 5(t), and h(r,t) is the funda- 
mental solution of the so-called non-Markovian forward drift equation 

u(t, t) = «o(r) - / K(t - s ) du ^ s ^ ds, r,f £ TZq , (6) 

JO OT 

where uq(t) = u(t, 0). 

When the kernel and the time-stretching functions are stated as 

K(t) = ^, g(t) = t a ^, 0<a<2, < /3 < 1 , (7) 

Equation (jl]) becomes 

P{X , t) - w + j^H jT r'l^ ^ - ^ ^ * , (8, 

that was originally introduced by A. Mura in his PhD Thesis [20J, and later 
discussed by him and collaborators in a number of papers [161 Ell [221 [23]. 

It is well-known that it exists a relationship between the solutions of a 
certain class of integral equations that are used to model anomalous diffu- 
sion and stochastic processes. In this respect, the density function P(x, t) 
which solves ([8]) is the marginal particle pdf, i.e., the one-point one-time den- 
sity function of particle dispersion, of the generalized grey Brownian motion 

(ggBm) gnuara. 

The ggBm is a special class of H-sssi processes of order H = a, or 
Hurst exponent H = a/2, where, according to a common terminology, H- 
sssi means H-self-similar-stationary- increments. The ggBm provides non- 
Markovian stochastic models for anomalous diffusion, of both slow type 
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< a < 1 and fast type 1 < a < 2. The ggBm includes some well-known 
processes, so that it defines an interesting general theoretical framework. In 
fact, the fractional Brownian motion appears for j3 = 1, the grey Brownian 
motion, in the sense of W. R. Schneider [301 [31], corresponds to the choice 
< a = f3 < 1, and finally the standard Brownian motion is recovered by 
setting a = j3 = 1. It is worth noting to remark that only in the particular 
case of the Brownian motion the stochastic process is Markovian. Moreover, 
the ggBm is not an ergodic process [22J. 

The integral in the non-Markovian kinetic equation ([5]) can be expressed 
in terms of an Erdelyi-Kober fractional integral. In fact, let fx, rj and 7 be 
fx > 0, rj > and 7 E 1Z, the Erdelyi-Kober fractional integral operator 
Z^ ,M , for a sufficiently well-behaved function (^(t), is defined as formula 
(1.1.17)] 

/.-J7O+7) r t 

r ^^ {t) = ~TW Jo rW1 {tV " rT ^ (r) d{jV) ' 

= JL-r^+i) / T'^-^f-T'rVMdr, (9) 



hence Equation (|8|) can be re-written as 



P(x,t) = P Q (x) + t c 



0,0 d^P 

a//3 Q x 2 



(10) 



The integral operator I™ was introduced by I. N. Sneddon (see for 
example [321 [331 El]) who studied its basic properties and emphasized its 
useful applications to the generalized axially symmetric potential theory 
(GASPT) and other physical problems (say in electrostatics, elasticity, etc). 
When rj = 1, one obtains the operators of fractional integration as originally 
introduced by H. Kober [9] and A. Erdelyi [1] and, when r/ = 2, those 
introduced by I. N. Sneddon [32j [331 El]- In the special case 7 = and 
77 = 1, the Erdelyi-Kober fractional integral operator ([9]) and the Riemann- 
Liouville fractional integral of order fx, here noted by J 1 *, are related by the 
formula 



<PV) = 44 / "(t - V(r) = r»J» <p(t) . (11) 



r(M) ^0 

The above remark about the relationship between the governing equation 
of the ggBm @ and the Erdelyi-Kober operator of Fractional Calculus ([9]) 
constitutes the aim of this Short Note. The possibility to re-write Equation 
([5]) as (|10|) was briefly noted by the Author in |25J. However, here this 
correspondence between the ggBm and the Erdelyi-Kober fractional integral 
operator is stressed and, since the ggBm serves as a stochastic model for the 
anomalous diffusion, this leads to define the family of diffusive processes 
governed by the ggBm as Erdelyi-Kober fractional diffusion. 
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In order to establish the diffusion- type equation corresponding to ([8|), 
the Erdelyi-Kober fractional differential operator is here introduced. Let 
n— 1 < fi < n, n E M, the Erdelyi-Kober fractional derivative is defined as 
[71 formula (1.5.19)] 



1 d 



(12) 



The Riemann-Liouville fractional derivative of order [a, m — 1 < \i < m, 

d m 

m £ M is defined as D^ipit) = ~^ l ^ m ^ </>(£); an d it emerges that the 
Erdelyi-Kober and the Riemann-Liouville fractional derivatives are related 
through the formula 



(13) 



A further important property of the Erdelyi-Kober fractional derivative is 
the reduction to the identity operator when \x = 0, i.e., 



(14) 



Recently, the notions of Erdelyi-Kober fractional integrals and deriva- 
tives have been further extended by Yu. Luchko [10J and Yu. Luchko & J. 
Trujillo p]. 

Equation ()10p in diffusive form is obtained by deriving in time both sides 
and it results 



dP 
~dt 



at c 



r0,/3 



d 2 p" 



dt \ a IP dx 2 



d 

a + t di 



a/P Q x 2 



(15) 



that can be recast as 



dP 
~dt 



« a-1 



09 



B d 

1+1 + ^- 

a at 



r 0,/9 



V/3 9 X 2 



(16) 



and finally, by using ([12 



(17) 



A diffusion-type equation for the ggBm was previously proposed [IB] but 
adopting, with an abuse of notation, the Riemann-Liouville fractional dif- 
ferential operator with a stretched time variable. Then, since the Erdelyi- 
Kober fractional differential operator is taken into account, Equation (|17p 
follows to be the correct formulation. 
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The Green function corresponding to (fTUl ED is [2(7| l2T| 122] [23] 

g ^ = \w^' 2 {^h) ' (18) 

where M v {z) is the M- Wright function, often referred to as Mainardi func- 
tion in the literature devoted to fractional diffusion [21 [26] , and it is defined 
as p3] 

oo {-z) n 

M v (z) = 2 , 

72=0 L V /J 

-j oc / \n— 1 

= - y / Txi T ( vn ) sin(vrw) , < v < 1 , (19) 

n=l v ' 

see Reference [21 El [TB] for a review. Here it is reminded the noteworthy 
composition, or subordination-type, formula [18J 



/-"A/„(!)=r* rM x ^)M e (L) with , A/. (20) 



where < z^,A,£ < 1 and £,t,r G 7£q- By using (|2U|) and the special 
case Mijiiz) = (1/\/tt) exp(— z 2 /4), Green function (fTB"j) can be expressed 

as [1511221123] 

SOM) = \^M» /2 (W\, (21) 



2 W2^/ 2 ' 

\ t -a/2 



I Z" 00 / \x \t~ a/z \ 

r°° l f x 2 1 

= X ^pM'^r^' (23) 

so that, under the view point of statistical mechanics, the ggBm, or the 
Erdelyi-Kober fractional diffusion, emerges to be the superposition of pro- 

1 f x 2 } 

cesses with stretched Gaussian density exp < — >, i.e. frac- 

tional Brownian motions, whose variance is (x 2 ) = 2rt a where r is a random 
coefficient distributed according to Mg(r). 

However, Equation (|23p can be further re-managed to exhibit a subor- 
dination type representation. In fact, after the change of variable i* = rt a , 
it follows that 



1 I X 2 1 1 . , ft 



e(XJ) = h wi-sj^W*' (24) 

which means that the random trajectory x = x(t) can be obtained as a 
subordination process by x = x(t) = y\t*{t)\, where t* = t*(t) is a positive 
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random variable that evolves in the natural time t and is referred to as 
operational time [HE]. The process i* = t*(t) is the directing process that 
realizes in the (t, i*)-plane whose pdf is t~ a M/3(t*t~ a ). Please note that the 
pdf of the directing process belongs to the same family of the Green function 
Q(x,t) and they differ for the parameter pair. The process y = y(t*) is the 
parent process that is a random trajectory in the (£*, y)-plane with Gaussian 
pdf evolving in the operational time i*. Geometrically, identifying the spatial 
coordinates y and x, the subordination structure x = x{t) = y[t*(t)] is 
obtained by concatenation. 

The marginal pdf of the non-Markovian diffusion process ggBm emerges 
to be related to the Mainardi function M v and it describes both slow and 
fast anomalous diffusion. In fact, the variance of Green function (|18p is 

/+oo 
x 2 Q(x, t) dx = (2/T(/3 + 1)) t a , then the resulting process turns 

out to be self-similar with Hurst exponent H = a/2 and the variance law is 
consistent with slow diffusion for < a < 1 and fast diffusion for 1 < a < 2. 
However it is worth noting to be remarked also that a linear variance growing 
is possible, but with non-Gaussian pdf, when f3 ^ a = 1, and a Gaussian 
pdf with non-linear variance growing when j3 = 1 and a/1. 

It is straightforward to note that, by using formula (fT3"j) . evolution equa- 
tion ([8]) reduces to the time- fractional diffusion if a = /? < 1, i.e., 

OP d 2 P 11 / Ixl \ 

^ =Z ^^' with Q ^ = 2m M ^m)^ (25) 

and variance (x 2 ) = (2/T(/3 + l))t^, and, by using formula (fT4|) . it reduces 
to the stretched Gaussian diffusion if a ^ 1 and (3 = 1, i.e., 

9P n 1 2 P . . 1 1 / \x\ 



11 x 2 \ 

ex P -77^ > ( 26 ) 



y/^TT t a / 2 \ At a 

and variance (x 2 ) = 2t a , and finally to the standard Gaussian diffusion if 
a = P = 1, i.e., (HJ) with ([2]) and variance (x 2 ) = 2t. The Green functions 
of these last two cases, i.e., (a ^ 1, j3 = 1) and (a = /3 = 1), follows by ([23]) 
noting that Afi(r) = S{r - 1). 

In general, even if the Green functions are interpreted as one-point pdf 
evolving in time, they cannot determine a unique (self-similar) stochastic 
process because this requires the determination of any multi-point pdf. But, 
what concerns the ggBm, since the increments are stationary, it emerges to 
be uniquely determined by its covariance structure |21|. l22j . Then, even if the 
ggBm is not Gaussian in general, it is a valuable example of a process defined 
only through its first and second moments, which indeed is a remarkable 
property of the Gaussian processes. Then the ggBm is a direct generalization 
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of the Gaussian processes and, in the same way, the Mainardi function M v is 
a generalization of the Gaussian function, and it emerges to be the marginal 
pdf of non-Markovian diffusion processes that describe both slow and fast 
anomalous diffusion. 

To conclude, in this Short Note it is highlighted the relationship between 
the Erdelyi-Kober fractional operators and the valuable family of stochastic 
processes generated by the ggBm, whose some remarkable properties are 
reported above, and the key role of the Mainardi function in this frame- 
work. In fact, the particle pdf of associated to the ggBm is the solution of 
a fractional integral equation (|10p . or analogously of a fractional diffusion 
equation (JTTJ) , in the Erdelyi-Kober sense and this solution is a Mainardi 
function. Since the governing equation of these processes is a fractional 
equation in the Erdelyi-Kober sense it is proposed to called this family of 
diffusive processes as Erdelyi-Kober fractional diffusion. 
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